Ground state and excitation dynamics in Ag doped helium clusters 
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We present a quantum Monte Carlo study of the structure and energetics of silver doped helium 
clusters AgHe n for n up to 100. Our simulations show the first solvation shell of the Ag atom to be 
composed by roughly 20 He atoms, and to possess a structured angular distribution. Moreover, the 
electronic 2 Pi/2 ^ 2 Si/2 and 2 P3/2 ^ 2 Si/2 electronic transitions of the embedded silver impurity 
have been studied as a function of the number of helium atoms. The computed spectra show a 
redshift for n < 15 and an increasing blueshift for larger clusters, a feature attributed to the effect 
of the second solvation shell of He atoms. For the largest cluster, the computed excitation spectrum 
is found in excellent agreement with the ones recorded in superfluid He clusters and bulk. No 
signature of the direct formation of proposed AgHe2 exciplex is present in the computed spectra of 
AgHeioo- 
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Superfluid 4 He clusters represent a gentle environment 
where high resolution spectroscopic studies of atoms, 
atomic clusters, and molecules at low temperature can 
be carried out |jl|. In such cold and fluid quantum sys- 
tems many perturbing effects due to the temperature and 
solid matrices are absent, allowing therefore for an eas- 
ier interpretation of the experimentally recorded spec- 
tra. Moreover, their superfluid behavior allows interest- 
ing quantum effects to take place and to be experimen- 
tally probed (for instance see Refs. Q, q])- 

Whereas the coupling of the rotational and vibrational 
motion of the molecules with the quantum motion of the 
solvent is permitted by the similarity between energy lev- 
els, the electronic structure of an atom is characterized 
by energy differences orders of magnitude larger than the 
ones needed to induce excitation in the atomic motion. 
Although this difference might seem to work in the di- 
rection of simplifying the physical description of the elec- 
tronic transition processes, many important details still 
wait to be clarified. As an example, the fluorescent D2 
emission line (i.e. the 2 Si/2 ^- 2 P3/2 radiative transition) 
of heavy single valence electron atoms dispersed in super- 
fluid helium is absent, while the Di line is sharp and only 
slightly shifted (1-2 nm) to the blue Q]. This is in con- 
trast with the large broadening and strong blueshift of 
the absorption lines. Moreover, some features of the LIF 
spectra of the dispersed Ag were interpreted as signature 
of the AgHe and AgHe2 exciplex formation || . 

The blueshift and broadening of the absorption lines 
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have been interpreted by means of a "bubble model". 
Here, the dispersed atom is enclosed in a spherical cavity 
due to the exchange repulsion of its valence electrons and 
the He ones. The liquid He around an atom is modeled 
by an isotropic sharp-edge density profile with no atomic 
internal structure. However, both the simple spherical 
bubble model Q and the one where quadrupolar distor- 
tions of the spherical cavity are allowed neither quan- 
titatively predict the absorption spectrum of Cs and Rb, 
nor allow to interpret the small splitting of the Rb D2 
line. Reasonably, the lack of any shell structure in the 
helium density profile, the absence of a full atomistic de- 
scription during the excitation process, and the physi- 
cally incomplete description of the bubble distortion by 
means of simple quadrupolar deformations may be held 
responsible for this undesired outcome [|| . 

In order to gain a better understanding of the excita- 
tion process and its dependency on the degree of " solva- 
tion" of the impurity, we feel a direct many-body simu- 
lation of the excitation spectra to be mandatory. This 
also allows to explore the change in the spectra upon the 
increase of the number of He atoms in the clusters, and, 
at the same time, to test the validity of our theoretical 
approach. 

With these goals in mind, we present a diffusion Monte 
Carlo study of the 2 P 3 / 2 $1/2 an d 2 Pi/2 ^~ 2 ^1/2 
absorption spectra of silver doped helium clusters. The 
Ag spectrum, both in bulk helium and in He clusters, 
has been deeply studied and well characterized 0, 0, [l(J 
showing that Ag is indeed solvated. Moreover, accurate 
interaction potentials between He and the excited 2 Si/2> 
2 Pi/2, an d 2 P3/2 states of Ag are available jy]. These 
potentials allowed to assign the broad band at 382 nm in 
the fluorescence spectrum to the AgHe2 exciplex |J. 

To tackle the atomic description needed to compute 
the excitation spectra, we believe the Monte Carlo meth- 
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ods are the best suited techniques. Since these meth- 
ods are well described in the literature |jl2|] , we sim- 
ply state that while variational Monte Carlo (VMC) al- 
lows one to optimize a trial wave function '^(R) and 
to successively compute any expectation values (O)vmc 
from it, diffusion Monte Carlo (DMC) corrects the re- 
maining deficiencies of the variational description pro- 
jecting out all the excited state components sampling 
/(R) = # (R)*r(R), or less commonly /(R) = *g(R). 

In atomic units, the Hamiltonian operator for our 
AgHe„ clusters reads as 

W = 4fe — + — WW (!) 
2 rru He m Ag J 

Here, we assume a pair potential of the form V(R) — 
J2i<j ^HeHe(^y) + J2i ^AgHcfVi) for the clusters with the 
silver atom in the 2 Si/2 electronic ground state. For 
VneHefrij) we employed the TTY potential jL3j, and for 
l / A g He(?'i) we fitted the 2 S AgHe potential by Jakubek 
and Takami We computed the energy of the AgHe 
dimer by means of a grid method ]l4j obtaining -4.021 
cm -1 : this value differs from their result (-4.000 cm -1 ) 
by only 0.021 cm" 1 . 

Our trial wave function has the common form 
<f T (R) = i>(nj) Elf 0( r i)i where no one-body part 
was used, and ip(r) = </>(r) = exp[— £§ — £J — E| — 
piT — Po ln(r)] + aexp[— b(r — tq) 2 ]. The parameters of 
the model wave function were fully optimized minimiz- 
ing the variance of the local energy for each cluster. The 
sampled distributions were used to compute exactly the 
energy using the mixed estimator 

{H)M = J /(R)dR (2) 

as well as the mixed and second order estimate (O)soe = 
2(C) m — (O)vmc of many other expectation values (e.g. 
the interparticle distribution functions) p2| . 

The resulting DMC energy values for the AgHe n clus- 
ters with n up to 100 are shown in Table | together with 
the differential quantity A(n) = — [E(n)— E(rn)]/(n— m), 
which can be interpreted as the evaporation energy of an 
He atom from the cluster. AgHe m is the largest cluster 
having m < n. 

From this data it can be noticed that A(rc) does not 
possess a monotonic behavior. Instead, the steady in- 
crease for n < 13 is followed by a rapid decrease in value 
before plateauing for n ~ 25. This behavior could be in- 
terpreted invoking different effects. For n < 13, a newly 
added helium feels the bare Ag interaction potential plus 
the interaction with the already present He atoms, that 
acts positively increasing the binding energy. Quantita- 
tively, we found the changes of A(n) versus n similar to 
the ones obtained for He„ pl|, He„H~ |jl6), and He„HF 
IjJ . Since this effect seems to be independent of the 
nature of the doping impurity, one may interpret it as 
a dynamical many-body effect of the interacting helium 
atoms. 



n 


E(n) 


A(n) 


Di 


D 2 


free Ag 






338.3 


328.1 
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-4.0212(9) 
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-8.2333(5) 


4.212(1) 344.5 


334.8 


3 


-12.598(2) 


4.365(2) 344.2 


334.6 


4 


-17.112(1) 


4.514(2) 






6 


-26.478(2) 


4.682(2) 342.6 


333.3 


8 


-36.259(4) 


4.890(2) 341.3 


332.3 


12 


-56.68(1) 




339.6 


329.8 


13 


-61.78(1) 


5.09(1) 






14 


-66.84(1) 


5.06(2) 






15 


-71.61(5) 


4.77(6) 


338.3 


328.5 


19 


-89.31(2) 








20 


-93.17(3) 


3.86(4) 


337.3 


327.6 


24 


-107.14(4) 


3.49(1) 






25 


-110.3(1) 


3.2(1) 






29 


-123.17(7) 








30 


-126.11(7) 


3.0(1) 


336.2 


326.6 


40 


-158.70(6) 


3.26(1) 


335.5 


325.7 


50 


-191.3(3) 


3.27(3) 


334.7 


324.9 


60 


-225.1(2) 


3.38(4) 


333.8 


323.9 


70 


-259.9(4) 


3.47(5) 


333.4 


323.6 


80 


-292.4(7) 


3.26(8) 


332.5 


322.7 


90 


-326.2(7) 


3.38(9) 


331.9 


322.0 


100 


-357.3(6) 


3.10(9) 


331.6 


321.7 



TABLE I: Total and evaporation energy (cm ), Di and D2 
absorption wavelengths (nm) for AgHe n clusters. 



Beyond AgHei3, the value of A(n) decreases indicat- 
ing the onset of a repulsive interaction. This could be 
attributed to an "excluded volume" effect, where each 
new He is strongly attracted by Ag in its first coordina- 
tion shell, but has to "find room" for itself forcing the 
other atoms to increase their local density, and rising 
their average kinetic energy. 

Finally, for clusters larger than AgHe25, the evapo- 
ration energy remains roughly constant around 3.1-3.5 
cm -1 indicating that a new He atom feels a quite differ- 
ent environment than for n < 25. 

Figure 1 represents the interparticle Ag-He probabil- 
ity density functions for clusters having n from 12 to 
100. These were normalized so that 47r r 2 p(r)dr = n, 
therefore representing the local density of He atoms 
around Ag. The functions for n < 12 overlap in shape 
with the AgHei2 one. As to AgHe3o, the presence of a 
broad shoulder at large r, that successively develops into 
a well defined peak, unambiguously indicates a second 
shell. More interestingly, the height of the first peak 
continuously rises until the second shell is completely 
filled, as indicated by the onset of another shoulder at 
large Ag-He distance for AgHeioo- Moreover, the density 
minimum between the first and second shell peaks also 
increases in height on going towards larger clusters, be- 
coming just 15% less than the second shell peak height. 
Both these evidences can be interpreted as a direct signa- 
ture of the "non rigidity" of the first He layer, as well as 
of an easy exchange process between the first and second 
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FIG. 1: He density distributions around Ag for n =12, 20, 30, 
40, 50, 60, 80, and 100. 



FIG. 3: Simulated absorption spectra for AgHe n clusters with 
n =12, 20, 40, 70, and 100. The vertical lines represent the 
free Ag spectrum. 
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FIG. 2: cos(HeAgHe) density distributions for n =2, 6, 8, 12, 
20, 30, 40, 50, 60, 80, and 100. Each distribution is shifted 
upwards with respect to the previous one by 0.1. 



shell |8|. 

As to the angular distributions, Figure 2 shows sev- 
eral cos(HeAgHe) distributions. The smaller clusters 
(n < 20) show a deep minimum for cos(HeAgHc) = 1 
and a smooth maximum located in the 0.6-0.8 range, 
both strong indications of a structured distribution of 
the He atoms in the first solvation shell. Here, the mini- 
mum indicates that two He atoms cannot overlap or sur- 
mount each other, the only two possible arrangements 
having cos(HeAgHe) = 1. Whereas the overlap is for- 
bidden by the repulsive part of the He-He potential, the 
possibility of an He atom to surmount another is hindered 
by the strength of the 2 S AgHe potential that forces 
the He motion in a limited radial region around Ag as 
shown by the AgHei2 radial distribution. Instead, the 
smooth maximum indicates the relative localization ef- 
fects due to the attractive interaction between He atoms. 



This effect is particularly evident for AgHe2, whose an- 
gular distribution function decreases on going towards 
cos(HeAgHe) = — 1. The position of the maximum shifts 
to larger cos(HeAgHe) values on going from n = 2 to 
n = 20, suggesting a progressively more structured pack- 
ing of the He atoms in the first shell for n — 20, and 
agreeing nicely with the aforementioned "excluded vol- 
ume" interpretation. The structured packing is also sup- 
ported by the shallow second peak located around 0.1 in 
the AgHe2o cosine distribution. Both the minimum and 
the maximum are "smeared out" by adding He atoms to 
AgHe2o, a clear indication that the second shell is less 
structured and more fluid than the first one. 

As to the absorption spectrum of the embedded Ag 
atom, we computed this observable using the semiclas- 
sical approach proposed by Lax |L9]], and adapted to 
the quantum Monte Carlo framework by Cheng and 
Whaley ^0|. In this method, the spectral lines are 
computed collecting the distribution of the difference 
V exc (R) — V gs (TV) or, more accurately, of the quantity 
V eX c{R) + J2i<j VkeHefaj) - E over the sampled /(R) 
fjlf . Here, V^ S (R) (V exc (R.)) is the interaction poten- 
tial between the ground (excited) state Ag atom with 
the sourronding He atoms, while Eq is the DMC ground 
state energy. The three V exc {R) PES for a given clus- 
ter configuration are obtained from the AgHe 2 Hi/2, 
2 n 3 / 2 , and 2 E interaction potentials using the Diatomic- 
in-Molecules approach. The details are well described by 
Nakayama and Yamashita for the Li, Na, and K cases 
II- 

The spectra obtained collecting V exc CR,) — V firs (R) dur- 
ing the simulations are shown in Table 0, and in Figure 3 
for several representative clusters. The same quantities 
obtained by collecting V exc (R) + J2i<j VkeHefaj) - E 
are blueshifted by less than 1 nm. The computed spec- 
tra clearly show the two separated bands deriving from 
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the excitation of Ag into 2 P\/2 and 2 i"3/2 states, the sec- 
ond one also displaying the classical short wavelength 
shoulder typical of the D2 line of heavy alkali atoms in 
superfluid helium |i[ Ffl]. For our largest cluster, the Di 
and D2 lines have maxima located at 331.6 and 321.7 
nm, and a FWHM of 4.3 and 9.8 nm, respectively. These 
results are in accurate agreement with the experimental 
wavelengths 332.8 and 322.5 nm, and FWHM 4.0 and 8.5 
nm @. 

From the spectra shown in Fig. 3, it clearly appears 
that the broadening of the absorption bands increases 
on going towards larger clusters. This evidence indicates 
that the Ag electronic degrees of freedom are coupled 
with the motion of an increasing number of He atoms, 
and not only with those located in the first shell. More 
interestingly, whereas all the clusters with n < 15 show 
a redshift with respect to the free Ag lines, the ones with 
n > 19 display a blueshift strongly dependent on the 
number of He atoms. Here, the redshift for n < 15 in- 
dicates that the clusters possess an internal distribution 
such that a vertical transition brings them in a region 
of the excited state potential where the complexes can 
form a bound state. This may give the possibility of 
producing AgHe„ (n =1-15) exciplexes starting from the 
corresponding clusters, and to experimentally study their 
spectrum and decaying dynamics. Conversely, the larger 
clusters are vertically excited to repulsive regions of the 
potential energy surface (PES), therefore preventing the 
direct formation of larger exciplexes. 



The blueshift for n > 19, at variance with basic solva- 
tion concepts, indicates a large effect of the second shell 
filling on the absorption wavelengths. This is confirmed 
by the computational evidence that the excitation spec- 
trum of AgHeioo, that shows the onset of a third shell, 
closely agrees with the one of AgHe 90 (see Table @). 

In order to rationalize this observation, as well as the 
monotonic blueshift of the absorption bands upon in- 
creasing of n, one must notice that the portion of the 
AgHe pair distribution located in the 10-13 bohr range 
overlaps with the tail of the repulsive excited AgHe 2 S 
potential. As a consequence, this zone of the pair density 
introduces some positive contribution to the diagonal el- 
ements of the matrix whose eigenvalues define the three 
electronic excited PES of the complexes. Since the mag- 
nitude of these contributions is dependent on the local He 
density via the sum V2 E (r) over the He atoms falling 
in that range, there is a net increase of the values of 
the diagonal elements upon increasing of the size of the 
cluster. This fact reflects itself in a positive shift of the 
eigenvalues, and hence in the blueshift of the computed 
spectra. 
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